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ABSTRACT 

Individual dark matter halos in cosmological simulations vary widely in their detailed structural 
properties such as shape, rotation, substructure and degree of internal relaxation. Recent non- 
parametric (principal component) analyses suggest that a few principal components explain a large 
fraction of the scatter in halo properties. The main principal component is closely linked with con- 
centration, which in turn is known to be related to the mass accretion history of the halo. Here we 
examine more generally the connection between mass accretion history and structural parameters. 
The space of mass accretion histories has principal components of its own. We find that the strongest 
two can be interpreted as the overall age of the halo and the acceleration or deceleration of growth at 
late times. These two components only account for ^ 70% of the scatter in mass accretions histories 
however, due to the stochastic effect of major mergers. Relating structural parameters to formation 
history, we find that concentration correlates strongly with the early history of the halo, while relax- 
ation correlates with the late history. We examine the inferences about formation history that can 
be drawn by splitting haloes into subsamples, based on observable properties such as concentration 
and shape at some final time. This approach suggests interesting possibilities, such as the possibility 
of defining young and old samples of galaxy clusters in a rigorous, quantitative way, or testing the 
dynamical assumptions of galaxy formation models empirically. 

Subject headings: cosmology: dark matter - large-scale structure of universe - galaxies: clusters: 
general - galaxies: groups: general 



1. INTRODUCTION 

Dark matter halos provide the framework for visible 
structure in the universe over a span of eight decades in 
mass, from the scale of rich galaxy clusters down to the 
scale of individual dwarf galaxies. Mass is normally as- 
sumed to be the main determinant of a halo's baryonic 
contents, and analytic models such as the halo occupa- 
tion distribution (HOD - Peacock & Smith 2000; Sel- 
jak 2000; Ma & Fry 2000; White et al. 2001; Berlind & 
Weinberg 2002; Cooray & Sheth 2002) make this assump- 
tion explicitly. Yet the halos that form in cosmological 
simulations vary greatly in shape, concentration, spin, 
substructure and other structural properties. As simula- 
tions of increasing size and resolution provide a more and 
more detailed picture of halo properties, and as observa- 
tional techniques including weak and/or strong gravita- 
tional lensing. X-ray and Sunyaev-Zel'dovich measure- 
ments (e.g. Corless et al. 2009; Oguri et al. 2010; Sereno 
& Zitrin 2011; Morandi et al. 2011) reach a precision 
where they can determine structural properties such as 
shape and concentration reliably, it is important to un- 
derstand how the structural features of a halo are inter- 
related, and what they can tell us specifically about its 
formation and evolution. 

Recently two groups, Skibba & Maccio (2011) and 
Jeeson-Daniel et al. (2011) (Sll and Jll hereafter), have 
taken the important step of performing non-parametric 
principal component analyses of halo properties. Prin- 
cipal component analysis (PCA) searches for simplifying 
trends in a complex data set by finding the axes in a 
multi-dimensional parameter space that account for the 
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largest fraction of the scatter. In the simplest case, it can 
uncover linear correlations in the data (e.g. fundamental 
lines or planes) and reveal hidden patterns or simplifica- 
tions. The results of Sll and Jll agree on some basic 
aspects of halo structural properties. Overall, the scatter 
in halo properties spans a fairly high-dimensional space, 
with 4 principal components required to explain about 
70% of the scatter. Nonetheless, a few strong compo- 
nents emerge. 

In terms of structural parameters, the first (strongest) 
principal component is best traced by concentration. 
The strong correlation with concentration suggests that 
this first principal component is linked to the overall age 
of the halo, and the strength of the correlation with zq.s 
confirms this. The origin of halo concentration has been 
studied extensively since it is a crucially important fac- 
tor in many calculations, including strong lensing (e.g. 
Broadhurst et al. 2005, 2008) and dark matter annihila- 
tion (e.g. Taylor & Silk 2003). Several analytic models 
have been developed over the years to explain concen- 
tration in terms of formation history (e.g. Bullock et al. 
2001; Wechsler et al. 2002; Zhao et al. 2009), so the link 
between formation history and this particular structural 
property is fairly well understood. Older systems, that is 
systems that had already assembled most of their mass 
into one or a few progenitors at early times, are more 
concentrated on average, although the exact details of 
the connection between age and concentration vary from 
one analytic model to another. 

If the first principal component of halo structural prop- 
erties is thus linked to age, we might naturally ask what 
the others correspond to. Are shape, spin or relaxation 
also related to the formation history, and if so how? To 
put this question in a quantitative framework, we first 
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have to decide how to describe the "forniation history" 
itself - what should we take this to raean exactly, given 
the complex set of merger and accretion events through 
which halos form? 

We can take the mass accretion history (MAH) as a 
starting point, defining this as the function A4{z) = 
M{z)/M{0) which describes the mass of the main pro- 
genitor of a halo as a function of redshift, normalized to 
the value at 2: = (van den Bosch 2002). Since M.{z) 
is a continuous function of a real variable, it contains an 
arbitrarily large amount of information about the history 
of a halo; equivalently, describing the MAH fully means 
specifying values of Ai{z) at an infinite set of redshifts. 
To characterize the MAH in simpler terms, we can turn 
once again to principal component analysis, approximat- 
ing each MAH as a vector of values M{zi) corresponding 
to the MAH evaluated at a finite, fixed set of redshifts zi. 
PC A of these vectors can then tell us whether mass accre- 
tion histories are well described by a set of basis functions 
characterzed by a single parameter, as suggested by van 
den Bosch (2002) and Wechsler et al. (2002), or whether 
they require two or more variables to explain their diver- 
sity, as suggested by Tasitsiomi et al. (2004) or McBride 
et al. (2009). 

Assuming a few principal components capture the 
main features of a halo's MAH, this will allow us to 
study correlations between structure and history in a 
well-defined and quantitative way. We note that this 
is only a first step towards understanding structure in 
terms of growth history; the MAH does not contain all 
the information about a halo's past by any means, and an 
alternate approach is to study the origin of a particular 
property such as shape or relaxation in detail, consider- 
ing the physical processes involved (e.g. Vera-Ciro et al. 
2011; Power et al. 2011). 

In this paper, we generate halos in a set of cold dark 
matter (CDM) simulations covering three different mass 
scales. We record the MAH of each well-resolved halo 
and decompose the space of MAHs, taken as vectors of 
values M{zi)^ into principal components. These compo- 
nents provide a non-parametric description of the MAH 
and clarify its basic properties. We also record the struc- 
tural properties of each halo at the present-day. We 
study the correlations between these properties them- 
selves (reproducing the trends found by Sll and Jll), 
and between structural properties and the principal com- 
ponents of the merger history. Finally, we discuss an ap- 
plication of our results, showing how samples of halos se- 
lected by concentration or shape will have systematically 
different formation histories. In this way, the observable 
properties of groups or clusters of galaxies can be used 
to infer their (unobservable) formation history. 

The outline of the paper is as follows. Section 2 de- 
scribes our n-body simulations and group finding, and 
how the structural properties of halos are defined and 
measured. Section 3 analyses the vector space of MAHs, 
decomposing them into principal components. In Section 
4, we then analyze the correlations between structural 
components, and derive principal components in this sec- 
ond vector space. In Section 5 we discuss applications of 
this work. We summarize our results in Section 6. 

Throughout the paper, we consider a WMAP7 cos- 
mology with parameters VLm = 0.27, I^a = 0.73, = 



TABLE 1 
Simulation Parameters 

Name Trip -^^box -^^start Ccom Cmax 

IO^QMq/i-i Mpc h-^ kpc h-^ kpc h'^ 

Sim60 0.01206 60 8L1 50 5 
Siml20 0.09648 120 67.3 110 11 
Sim240 0.77184 240 55.2 240 23 

(in our CDM-only simulations), = 72 kms~-^Mpc~"^, 
Us = 0.96, and erg = 0.80. 

2. MEASURING HALO PROPERTIES 
2.1. Numerical Simulations 

Our simulations were performed using the massively 
parallel n-body code Gadget2 (Springel 2005), with cos- 
mological initial conditions generated by initial condition 
code Graf ic2 (Bertschinger 2001). Three different box- 
sizes were simulated with Np = 512^ particles each, pro- 
ducing halo catalogs that vary in mass by a factor of 8 
from one simulation to the next. 

The softening length used in the simulations was 
e = min(ecom7 ^max/^i), where ecom is a comoving soft- 
ening length and 6niax 

is a (maximum) physical soften- 
ing length. The softening lengths and other simulation 
parameters are listed in Table 1. 

In each simulation, we output 95 snapshots spaced 
logarithmically in cosmic scale factor a between z = 8 
and z = 0. To detect halos in each snapshot we 
use the Friends-of-Friends (hereafter FOF) algorithm 
(Davis et al. 1985) with a linking length of b = 0.2 

times the mean inter-particle separation LboxN"^^^. 
We also use the SMOOTH algorithm (Stadel 1995) 
to compute the mean local density around every par- 
ticle. (These tools can be found at http://www- 
hpcc.astro.washington.edu/tools/tools.html.) From the 
initial FOF catalog we select all halos with 1000 particles 
or more, reducing the number of usable FOF groups to 
1560, 1580, and 1532 respectively for simulations Sim60, 
Siml20, and Sim240. This gives us a total sample of 4672 
well-resolved halos. 

In a sample defined solely by FOF linking, many 
groups will be multi-component systems caught in the 
act of merging. This can introduce significant scatter in 
the measured structural properties of halos. We can flag 
unrelaxed systems by using the minimum value of Xnfw 
(see Eqn. (3) below) as a goodness- of- fit indicator. Al- 
ternately, the asymmetry of the mass distribution can 
also be used to distinguish relaxed and unrelaxed sys- 
tems. The center of each halo is taken to be the position 
of the particle with the highest local mean density within 
the FOF group, Xp, as centering on this point generally 
appears to produce smoother fits to the density profile 
than if we use the center of mass of the group, xcoU- 
The offset between Xp and xqoM relative to the size of 
the halo is then an measure of asymmetry or relaxed- 
ness. Defining Xoff = \xp — ^coMl/^vir, we select from 
our full sample a sub-sample of 3290 relaxed halos for 
which Xoff < 0.07 and Xnfw ^ O-^- These correspond 
to the criteria used by Maccio et al. (2007) and Sll (our 
Xnfw their Prms)- We note that they are fairly inclu- 
sive, defining 70% of our sample as relaxed; Power et al. 
(2011) suggest using the stricter criterion Xofr < 0.04, for 
instance, to select relaxed systems. In what follows, we 
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will consider results both for the full sample and for the 
relaxed sub-sample. 

2.2. Density Profiles 

To determine the spherically-averaged density profile, 
we bin the halo particles radially in equal-sized bins, 
starting from the center of the halo (i.e. Xp) and pro- 
ceeding outwards until we reach the virial radius Tvir, 
defined by 

Mvir = (47rr^ir/3)Avir,cPc , (1) 

where Mvir is the virial mass of the halo, taken to be the 
FOF mass, and Avir,c = 200. We note that a fraction of 
the outer particles linked by FOF may lie outside Tvir, 
particularly at low redshift. 

We fit a Navarro-Frenk- White (Navarro et al. 1996, 
NFW hereafter) density profile to each halo over the ra- 
dial range [O.Olrvir, ^vir]- We do not include the inner- 
most 1% of the virial volume in the fit, since resolution 
effects and softening may systematically alter the density 
profile there. The NFW profile is 



Pnfw(^) 



(r/r,)(l + r/r,)2 



(2) 



where the scale radius (the radius at which the loga- 
rithmic density slope is —2) and ps is the characteristic 
density. Best fit values of these parameters are obtained 
by minimizing the quantity 

Xnfw = ^ \^^P(^^) - InpNFV^(ri)]^ , (3) 

i^bins .^^ 

where Nbins is the number of bins used in the fit. This 
minimization technique is logarithmic rather than linear 
in order to give similar weight to the fitting near the 
centre and at the outer edges of the halo (Jing 2000; 
Maccio et al. 2007). The minimization code uses the 
Levenberg-Marquardt method to find the best fit values 
of Ts and ps. As mentioned previously, we find that using 
the particle with the highest local density as the center of 
the halo produces good NFW fits to the density profile, 
while fits using the centre of mass have larger residuals. 
Bin size also affects our results somewhat; after some 
experimentation we chose a bin size of 15 particles to 
measure the profile and determine halo concentration. 

2.3. Concentrations 

Given an NFW density profile fit, the concentration 
parameter is defined as c = Tvir/^s- Here rvir is the virial 
radius, that is the radius around the halo center that 
encloses mass Mvir, as defined in Eqn. (1). Physically, 
concentration is a measure of the size or density of the 
central core relative to the size or density of the whole 
halo. 

We could determine the concentration of a halo by 
measuring Tyir and r^ separately and taking the ratio 
of the two. In practice, a slightly different technique ap- 
pears to give more robust results. Integrating the NFW 
profile with respect to radius out to rvir, the enclosed 
mass is given by 



Mvir = ^TTpsrl 



ln(l + c) - 



1 



(4) 



Comparing the definition of the virial mass (Eqn. (1)) 
with Eqn. (4), we can derive the non- linear relation 



/(c) 



1 



ln(l + c) 



1+c 



r,cPc 



3ps 



(5) 



Solving this equation then gives us the concentration c, 
as in (Klypin et al. 2002; Gavazzi et al. 2003). We wih 
use this second method throughout the rest of the pa- 
per, although we find the two methods for computing 
concentration are generally consistent with one-another 
to within 2-3%. Typical errors in concentration are large 
(20-30%) for our smallest halos, but drop below 10% for 
halos with more than 3,000 particles. 

2.4. Mass Accretion Histories and Ages 

To determine the MAH for a given halo, we construct 
its merger tree, working backwards from the final snap- 
shot at ^ = 0. The merger tree links the final halo to its 
progenitors at each previous redshift step. To determine 
whether a halo Hiater found at a later time is related to 
an earlier halo Heariier, we require that more than half 
of the particles in 'Heariier be contained in 'Hiater- If this 
condition is satisfied then we say that halo 'Heariier is a 
parent halo (or progenitor) of the child halo Hiater- This 
method of defining parent halos restricts all earlier halos 
to have at most 1 child halo in the consecutive snapshot. 
The method does not restrict the number of parent ha- 
los a child halo may have, so in order to construct the 
MAH, we define the main parent to be the one with the 
largest contribution of particles to the child halo. This 
then produces a single sequence of parent halos in the 
merger tree, so we can assign a well-defined MAH to the 
halo in the final snapshot. 

We fit our final MAHs with the 2-parameter model 
proposed by Tasitsiomi et al. (2004) and McBride et al. 
(2009) 



M{z) = M{z)/M{0) = (1 -h zfe- 



(6) 



where /3 and 7 are free parameters to be fit, and M(0) is 
the mass of the halo at 2: = 0. This form is a generaliza- 
tion of the earlier 1 -parameter fit proposed by Wechsler 
et al. (2002). The physical interpretation of the param- 
eters f3 and 7 can be seen by noting that 



dz 



lnM{z) = -(1 + zfe-^' [I3{1 + z)'^ - 7] (7) 
= 'y-l3 for^-0 (8) 



Thus (7 — /3) is the logarithmic growth rate of the halo 
in the limit z ^ (approaching from above). 

We perform the model fits to each MAH by minimizing 
the quantity 



Xm 



1 



N« 



snaps 

E 

i=l 



M(0) 



■M{zi) 



(9) 



where Ngnaps is the number of snapshots for which the 
halo has existed in our sample. 

2.5. Shape and Triaxiality 
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In order to describe the shape of a halo, one can diag- 
onalize its moment of inertia tensor 

N 

^ij = ^rnn {\\rn\\^Sij - ri^nrj,n) (10) 

n=l 

where r^^^ is the i*^ Cartesian component (relative to 
the halo center) of the n}^ particle of the halo. This will 

produce eigenvalues Xf that are related to the relative 
axis lengths {a, 6, c} of the a halo (assuming an ellipsoidal 
mass distribution) by 

Af = ^Mhalo(6' + c2) (11) 

xT = ^Mhaio (a' + c2) (12) 
5 

xf = ^Mhaio (a' + 6') . (13) 

With simple linear combinations of Eqn. (11) - (13), one 
can then obtain the relative axis lengths. 

A more direct - but equivalent - method of obtain- 
ing the relative axis lengths {a, 6, c} is to diagonalize the 
tensor describing the second moment of the mass distri- 
bution 

N 

^ij = ^ rnnVi^nrj^n (14) 

n=l 

The 3 eigenvalues Xf^ from the diagonalization of the 
matrix in Eqn. (14) will be the squares of the relative axis 
ratios, so the relative lengths of the principal axes can 

thus be determined as {a, 6, c} = {|A^ I^J^f I^J-^f 1^}- 
By convention, throughout the rest of the paper we will 
we re-label the axis lengths as necessary so that a > 
b > c. We note there are several other possible ways to 
measure shape - see (Zemp et al. 2011) for a discussion 
of the relative advantages of different techniques. 

As well as taking ratios of the individual axis lengths, 
we will also consider the triaxiality index T = {o? — 
6^)/(a^ — c^). Oblate spheroids have T = 0, while purely 
prolate halos have T = 1. Halos tend to have a skewed 
distribution of T values, so we also consider the elonga- 
tion parameter E = {b'^ -\- c^)/(2a^), whose distribution 
is more Gaussian. 

2.6. Summary of measured halo properties 

In studying correlations between halo properties, we 
will consider a large set of properties so as to overlap 
with the two recent PC A studies (Sll and Jll). Clearly, 
not all the parameters listed below are independent. For 
example, we will assess several different definitions of 
shape, but they all ultimately depend on the two axis 
ratios c/a and b/a. Our main goal will be to reproduce 
the patterns seen in Sll and Jll, and then relate them 
to aspects of the MAH. Any redundant parameters will 
appear as strong correlations in our analysis, and can be 
combined or ignored accordingly. 

1. Virial Mass Mvir 

This is the total mass of halo within a virial radius 
^vir = '^200, c of the halo centre, which is taken to 
be the particle of highest local density within the 



linked group found by FOF. The virial radius is 
defined such that the halo has an average density 
of 200pc. 

2. Concentration C200 

The concentration c = Vyir/f^s is determined using 
Eqn. (5), as explained above. 

3. Formation Redshift Zx for x G (0, 1) 

This is the redshift at which the MAH reaches 
M{z) = M{z)/M{0) = X. The most common ex- 
ample is Z0.5, at which M(zo.5) = 0.5M(0). 

4. Mass Fraction History {M/Mq)^ for z > 

The mass fraction history is a complementary mea- 
sure of age, equal to M{z) for a specified value of 
z. 

5. Relaxedness XofF 

This is the distance between the center of mass of 
the halo and the particle of highest local density, 
divided by the virial radius of the halo Tvir- It 
is expected that relaxed halo without any recent 
mergers will have a small value for Xofr- Through- 
out this paper, we will use Xofr < 0.07, Xnfw ^ ^-^ 
as a relaxedness condition. 

6. Triaxiality T 

Defined as T = (a^ - 6^) /{o? - c^) where a> b> c 
are the lengths of the principal axes of the halo. T 
measures the prolateness or oblateness of a halo. 
Spherical and sausage-shaped halos have T 1, 
while disk-shaped halos have T ?^ 0. 

7. Elongation E 

This is defined as E = [{b/af + {c/af]/2. The 
distribution of values of E is more Gaussian than 
the distribution of triaxiality values. 

8. Sphericity c/a 

This is a measure of the sphericity of the halo, di- 
viding the shortest axis length c by the major axis 
length a. 

9. Spin A 

To provide a measure of rotation, we use a dimen- 
sionless spin parameter A = Jvir/(2GM^-^rvir)^^^ 
where Jvir, Mvir, and 

Tvir are the total angular 
momentum, virial mass, and virial radius, respec- 
tively. Discussion of a different spin parameter def- 
inition can be found in Maccio et al. (2007). 

10. Environment D^j for n eZ and / > 

The dimensionless environment parameter Dnj is 
defined as distance of the n*^ nearest halo that has 
a virial mass greater than /Mvir, divided by its 
virial radius Tvir- This definition is physically mo- 
tivated because it scales as the tidal force to the 
-1/3 power (Haas et al. 2011). 

3. ANALYSIS OF THE MASS ACCRETION HISTORIES 

3.1. The (13, -f) Fit 

For each of our well-resolved halos, we fit the func- 
tional form in Eqn. (6) to the MAH to determine values 
of p and 7. Fig. 1 shows the resulting distribution of pa- 
rameter values. The colour coding indicates the value of 
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Fig. 1. — Fitted values of 7 versus (3 for 4672 halos with well-defined MAHs. The color scale indicates the value of relaxedness parameter 
Xoff? with red points being more relaxed (left-hand panel), or the formation redshift log]^Q(2;o.5), with red points being older (right-hand 
panel). The large points on the right-hand panel indicate the values of /3 and 7 used to fit the curves in Fig. 2 (triangles - PCI; squares - 
PC2). 



the relaxedness parameter Xoff (left-hand panel), or the 
formation redshift as log (^0.5) (right-hand panel). 

While this functional form of Eqn. (6) produces a rea- 
sonable fit to most of our MAHs, we can see immediately 
that p and 7 are not necessarily a natural choice of vari- 
ables, as there is a strong correlation between them for 
most MAHs. Furthermore, these parameters are degen- 
erate in many cases, with different combinations of values 
providing almost identical fits to the MAH, as discussed 
in (Taylor 2011). A more interesting variable is the quan- 
tity (7 — /3) , which corresponds to the logarithmic growth 
rate as z ^ 0. Given the simplicity of the (/3,7) fit, this 
quantity also fixes the overall age of the system, as can be 
seen clearly from the color scale in the right-hand panel 
of Fig. 1. We also indicate on the right-hand panel of 
Fig. 1 the values of P and 7 used to fit the curves in 
Fig. 2 (see below). 

3.2. Principal Component Analysis of the MAHs 

Principal Component Analysis (PCA) is a non- 
parametric technique for decomposing a set of Np possi- 
bly correlated variables into an equal or smaller number 
of independent variables. Because the final number of 
variables may be smaller than the initial number, PCA 
has the ability to reduce the dimensionality of a data set 
and uncover hidden patterns. 

Technically, PCA involves determining the eigenvec- 
tors (principal axes PCi) of the covariance matrix and 
their respective eigenvalues (principal components A^). 
Using the convention that principal components are ar- 
ranged by magnitude such that Ai > A2 > ... > Aatp, 
the principal components A^ are the relative fraction of 
the total variance in the direction of the principal axes 
PCi. Thus PCA provides us with principal axes PCi 
pointing in the direction of greatest variance in the data, 
PC2 pointing in the direction of second greatest variance 
orthogonal to PCi, etc. Because the individual proper- 
ties analyzed by PCA may have different dimensionality 
and/or magnitude, it is important to transform the ini- 
tial data such that the mean of each parameter is and 
the variance is 1, assuring that all the input fields will be 
treated equally. Standardizing the variables in this way 
will also scale the principal components to sum to Np, 



so the relative contribution of the principal component 
Xi to the variance is the fraction Xi/Np. 

We can use PCA to derive a more fundamental decom- 
position of the MAHs, taking as the input data vectors of 
values M{zi) for a specified set of Zi. The resulting prin- 
cipal axes will then capture the basic shape variations 
of the MAHs. For this analysis we consider a subset of 
MAHs for 817 halos extracted from Siml20 that are in 
the well-resolved sample (1000 particles or more, corre- 
sponding to masses of 10^^ or more) at z = 0, and 
have parents detected above the basic FOF limit of 200 
particles or more (corresponding to masses of 2 x 10^ ^Mq 
or more) for at least 70 outputs (i.e. going back to 2: ~ 4). 

Fig. 2 shows the mean MAH for these systems (cen- 
tral black curve), plotted in terms of the scale factor 
a. The first principal component (MAH-PCi), which 
accounts for 51% of the total variance in the MAHs, rep- 
resents the most important variations away from this av- 
erage. We can determine the nature of MAH-PCi by 
adding or subtracting from the mean MAH the first prin- 
cipal component vector, normalized to 1 or 2 times the 
r.m.s. scatter along the MAH-PCi axis. The resulting 
curves, illustrating ±1-2 a deviations along the MAH- 
PCi axis, are indicated by the upper (red) and lower 
(blue) curves in the left-hand panel of Fig. 2. Clearly, 
MAH-PCi corresponds to the overall age of the halo, 
with the upper (red) curves indicating older systems and 
the lower (blue) curves indicating younger systems. Both 
the mean and the ±1-2 a deviations are all well fit by 
the McBride model (gray curves). The best fit values of 
P and 7 are plotted as triangles on the right-hand panel 
of Fig. 1. Roughly speaking, they are spread out in the 
(7 ± /3) = constant direction, but the precise positions 
appear complicated in this parameter space. This may 
be the result of the degeneracies in the (/3, 7) fit discussed 
earlier. 

The second principal component, MAH-PC2, ac- 
counts for a further 17% of the variance. The right-hand 
panel of Fig. 2 shows the effect of ±1-2 a deviations 
away from the mean along the MAH-PC2 axis. This 
variation is more complicated than that of MAH-PCi, 
but roughly speaking it represents acceleration or de- 
celeration in the MAH. The (/3, 7) model is no longer 
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Fig. 2. — The mean MAH for 817 halos from Siml20 (black solid line), and its variation as one moves zbl and ±2 standard deviations 
along the first (left-hand panel) and second (right-hand panel) principal component axes. In the left-hand panel, the upper (red) curves 
represent older halos, and the lower (blue) curves represent younger halos. In the right hand panel, they represent halos whose accretion 
rate has accelerated or decelerated close to z = 0, respectively. The gray lines indicate McBride fits to the curves. 

a particularly good fit to these curves, particularly the 
ones that accelerate rapidly at late times. The relation- 
ship between the best fit values of the P and 7 is much 
simpler in this case, as indicated by the squares on the 
right-hand panel of Fig. 1, but given the poor fit to the 
curves it is not clear whether this is significant. 

We have examined the other principal components of 
the MAH. Generally speaking, they appear to be Fourier- 
like decompositions of the curve with increasing numbers 
of oscillations about the mean MAH. Given that individ- 
ually they account for relatively little of the scatter (9% 
for MAH-PC3, 6% for MAH-PC4, < 5% for subse- 
quent PCs), we will not discuss them any further. We 
note however that the total power in PCs 3 or higher is al- 
most 30%, showing that MAHs are not particularly well 
fit by smooth curves, presumably because of the large 
jumps produced by major mergers. 

It is also worth noting that the distribution of PCA 
vector components for individual MAHs is very far from 
a Gaussian scatter about the mean. Fig. 3 shows the dis- 
tribution of the MAH-PCi coordinate values for the 817 
MAHs considered in Fig. 2. Most MAHs have either high 
or low values of the MAH-PCi component, suggesting 
that the mean MAH is not necessarily very indicative of 
individual MAHs, which tend to deviate strongly from 
the mean. 




3.3. Relation to Excursion Set Theory 

Finally, we can ask to what extent the basic shape 
of the mean MAH, or the variations of MAH-PCi 
around it, reflect the statistics of halo growth predicted 
by excursion-set (or extended Press- Schechter - EPS) 
theory. In excursion set theory, the growth of a halo 
can be considered a random walk in the space of vari- 
ables cr^(M) and uj{z), where <j^(M) is the variance of 
the primordial density field smoothed on a mass scale M, 
and uj{z) = Sc{z) is the height of the barrier for spher- 
ical collapse at redshift z (Lacey & Cole 1993). The 
distribution of jumps Aa^ over a redshift step corre- 
sponding to Auj is a simply a Gaussian in (Acj/a/Ao^) 
(cf. Lacey & Cole 1993, Eqn. 2.15). This suggests that 
MAHs should look much simpler in the coordinate sys- 
tem (cj, Acr^). Fig. 4 shows the mean MAH and the 



-0.5 0.5 1 

Projection onto MAH-PC1 

Fig. 3. — The distribution of the MAH-PCi coordinate value 
for 817 weh-resolved halos from Siml20. 

variations along the MAH-PCi axis in (cj, Aa^) coor- 
dinates. Sure enough, the different curves from Fig. 2 
are now almost straight lines in these new coordinates. 
On closer examination, there is still a slight residual cur- 
vature to the lines. Tests with artificial MAHs suggest 
this is at least partly due to the difference between the 
variables Aa^ = a'^{Mj) —a'^{Mo), which is plotted here, 
and Acr^ = cr'^{Mj) - cr'^{Mi), which is assumed in EPS 
theory, for a section of the trajectory going from Mj to 
Mi > Mj . It may also refiect from corrections introduced 
to EPS statistics by ellipsoidal collapse (Sheth et al. 2001; 
Sheth & Tormen 2002). Overall, however, the transfor- 
mation from (cj, cr^) to (a, M) coordinates accounts for 
most of the characteristic curvature of MAHs in (a, M) 
space. This is useful, as it gives an approximate way 
of predicting MAH distributions as a function of back- 
ground cosmology. 

4. ANALYSIS OF STRUCTURAL PROPERTIES 

4.1. Correlations Between Structural Parameters 

The simplest way to search for relationships between 
the structural properties of halos is to measure the cor- 
relation of one property with another. We will use the 
Spearman Rank coefficient as a quantitative measure of 
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Fig. 4. — The mean MAH (black line) and its variation by ±lcr 
and zb2cr along the first principal axis (upper/blue and lower/red 
solid lines), plotted in terms of the EPS variables Acj and Acr^. 

the correlation strength. The Spearman Rank coeffi- 
cient, rsRc ^ [—1, 1], is a non-parametric measure of cor- 
relation which determines the strength of the relation- 
ship between two variables based on whether the data 
can be characterized by a monotonic function. Perfect 
Spearman correlation corresponds to rsRc = 1 (positive 
correlation) or rsRc = — 1 (negative correlation), while 
^SRc ^ indicates uncorrelated variables. 

The full correlation analysis of all 11 halo properties is 
shown in Appendix A. Broadly speaking, the strongest 
correlations are between age indicators such as the for- 
mation redshift logzo.5, and concentration c. Many pa- 
rameters, such as logzo.5 and log[A^(0.8)], are nearly de- 
generate. Some weaker but physically significant corre- 
lations exist between mass and age indicators, and be- 
tween shape and concentration. The former is easily un- 
derstood as a consequence of hierarchical structure for- 
mation (more massive halos have formed more recently). 
The latter appears to be a newly discovered pattern. 

Correlations between halo properties have already been 
analyzed by Jll. In order to compare directly with their 
results, we also show, in Appendix A, a table similar to 
theirs. We omit substructure because we lack the res- 
olution to detect it convincingly in many of our halos. 
Overall our results are almost identical to theirs; concen- 
tration and age are the most fundamental parameters, 
with mass and shape following behind. We find as they 
do that environment is not strongly correlated with any 
of the halo properties. The only result that differs signif- 
icantly is the correlation between the spin parameter A 
and concentration. It is not clear what the origin of this 
discrepancy is. 

4.2. Principal Component Analysis 

We perform PCA, as described in Section 3.2, on the 
11 halo properties described in Section 2.6. Our results 
are broadly consistent with the previous analyses of Sll 
and Jll, but we include them here for completeness. The 
strongest structural principal components (S-PCs) with 
their respective principal axes are listed in Table 2 for 
the full sample and Table 3 for the relaxed subsample. 
Numbers in bold indicate significant contributions to the 
principal component, taken here to be coordinate val- 
ues \vi\ > 0.3. We note that several of our parameters 



are almost degenerate, so our PC vectors have signifi- 
cant contributions from more parameters than in previ- 
ous analyses (Sll, Jll). We summarize the properties of 
the four strongest PCs below. 



1. S-PCi includes concentration and the various age 
indicators (zo.2, ^o.s? (^/^o)o.5) as its strongest 
contributors. Mass and elongation also contribute 
significantly to PCI for the relaxed sample; for 
the whole sample relaxedness Xoff is more impor- 
tant than mass. The overall weight of S-PCi and 
its composition are roughly in agreement with the 
PCA from Sll and Jll; we verify that concentra- 
tion, age, mass, and shape are seem to be related, 
fundamental parameters. 



2. S-PC2 contributes half to a third as much variance 
as S-PCi, and includes significant contributions 
from shape, age, and spin. In fact, the contribution 
from shape parameters (T, c/a) is stronger than 
in S-PCi. This suggests that some aspect of the 
MAH of dark matter halos helps determine their 
shape at 2; = 0. 



3. S-PC3 is then fairly distinct from the previous 
components, depending mainly on mass and spin. 
In the case of relaxed halos, environment also con- 
tributes at a lower but significant level to S-PC3. 



4. S-PC4 consists mainly of environment, showing 
that environment is a fairly independent param- 
eter, or conversely that structural properties are 
largely independent of environment. For the re- 
laxed subset there is some contribution from re- 
laxedness and triaxiality, showing these do corre- 
late with environment in relaxed systems. 



TABLE 2 

Principal components of the full sample 



Principal axis 


S-PCi 


S-PC2 


S-PC3 


S-PC4 


X^/Xl 


1.00 


0.34 


0.27 


0.22 


Contribution 


40% 


14% 


11% 


9% 


log Mvir 


0.29 


0.11 


0.53 


-0.08 


log C200 


-0.39 


0.04 


-0.04 


-0.13 


log 2:0.2 


-0.39 


0.29 


-0.10 


0.06 


log 2:0.5 


-0.39 


0.32 


0.01 


0.12 


log (M/Mo)o.5 


-0.37 


0.33 


0.01 


0.16 


log Xoff 


0.30 


-0.05 


-0.10 


0.26 


logT 


0.22 


0.42 


-0.21 


0.21 


logE 


-0.34 


-0.46 


0.21 


-0.12 


log c/a 


-0.27 


-0.41 


0.08 


0.05 


log A 


0.03 


-0.28 


-0.76 


0.01 


log Di^o.i 


0.07 


0.23 


-0.15 


-0.90 



8 



Wong & Taylor 



TABLE 3 

Principal components of the relaxed sample 



Principal axis 


S-PCi 


S-PC2 


S-PC3 


S-PC4 




1.00 


0.43 


0.29 


0.24 


Contribution 


36% 


16% 


10% 


9% 


log Mvir 


0.34 


0.11 


-0.46 


-0.08 


log C200 


-0.39 


0.10 


0.19 


-0.06 


log 2:0.2 


-0.41 


0.27 


0.06 


0.02 


log 


-0.40 


0.32 


-0.06 


0.08 


log (M/Mo)o.5 


-0.36 


0.30 


-0.11 


0.17 


log XoR 


0.17 


-0.16 


0.04 


0.36 


logT 


0.21 


0.36 


0.20 


0.43 


log E 


-0.33 


-0.49 


-0.20 


-0.15 


log c/a 


-0.28 


-0.43 


-0.13 


0.08 


log A 


0.01 


-0.31 


0.72 


0.19 


log 1^1,0.1 


0.10 


0.20 


0.35 


-0.76 



4.3. The origin of halo shapes 

PCA on the structural properties of halos in Section 
4.2 indicated correlations between formation history and 
shape, particularly in (structural) PC2 (cf. Tables 2 and 
3). There are various ways in which formation history 
could influence shape. One simple possibility is that the 
final shape of a halo is determined by the shape of the 
initial, uncollapsed matter from which it forms. This 
connection would not be apparent in the MAH, since it 
does not record the shape of the merging material, but 
only the rate at which it merges. We can easily test 
the hypothesis that final shape tracks initial shape by 
measuring the latter in early timesteps, including all the 
particles which eventually merge to form a halo. 

Fig. 5 shows how the elongation of a halo at 2: = cor- 
relates with the elongation of the particle distribution 
from which the halo formed, measured at a higher red- 
shift, for the full (left-hand pane) and relaxed (right-hand 
panel) samples. There is a strong correlation between ini- 
tial and final elongation looking back to recent redshifts, 
particularly for relaxed halos; at higher redshift or for less 
relaxed halos the correlation gradually weakens. Table. 4 
gives the Spearman rank coefficients iFor the correlations, 
showing that they are significant back to z = 1. The 
implication is that halo shape can trace the geometry of 
recent phases of the merger process. We note that Vera- 
Ciro et al. (2011) have recently completed a much more 
thorough study of the nature and mechanisms by which 
shape evolves as a halo grows. 

TABLE 4 

Correlation Between Final and Initial Elongation 



z 


ve (all halos) 


ve (relaxed halos) 


0.1 


0.630 


0.727 


0.5 


0.338 


0.444 


1.0 


0.300 


0.405 


2.4 


0.082 


0.128 



4.4. Structural Properties Versus MAH-PCs 

Given our earlier non-parameteric decomposition of 
the MAHs into their own set of principal components 
in Section 3.2, we can examine how structural prop- 
erties relate to the main features of the MAH. Fig. 6 



shows how Z0.5, concentration and relaxedness Xoff corre- 
late with MAH-PCi and MAH-PC2, for the fuh sam- 
ple (left-hand panels) and the relaxed sub-sample (right- 
hand panels). In the first row of panels, we see that the 
relationship between 2:0.5 and the two MAH-PCs is com- 
plex. MAH-PC2 is fairly cleanly ant i- correlated with 
Z0.5, in the sense that accelerating MAHs have lower val- 
ues of 2^0.5, but the relationship between MAH-PCi and 
2^0.5 more complicated. For large values of zq.s, the two 
parameters are reasonably ant i- correlated, but systems 
with low values of 2:0.5 show a large scatter in 2:0.5 at 
almost constant MAH-PCi ^ 0. These are probably 
systems with recent major mergers in their MAH, al- 
though it is interesting that many are present even in 
the relaxed subsample (right-hand panels) . This may in- 
dicate that the criteria we have used to define the relaxed 
sample (xofr < 0.07 and Xnfw < 0-^) ^^o inclusive. 
Power et al. (2011), for instance, use the stricter criterion 
Xoff < 0.04. 

Concentration (second row of panels) follows a simi- 
lar pattern. At high concentration, concentration and 
MAH-PCi are well correlated, while at low concentra- 
tion there is a large scatter at almost constant MAH- 
PCi. In this case, however, these are often unrelaxed 
systems; restricting ourselves to the relaxed sample, 
we see a cleaner correlation between concentration and 
MAH-PCi. Overall, the distribution suggests a simple 
relation between concentration and MAH-PCi, which 
is occasionally disturbed by mergers. A similar ana- 
lytic model relating concentration to the MAH has been 
proposed by Zhao et al. (2009); they distinguish phases 
of slow growth, where concentration gradually increases, 
from phases of rapid growth, when concentration is re- 
set to a low value. Thus concentration should be a good 
age indicator for old, relaxed systems, but a poor one for 
low-concentration and/or unrelaxed systems. 

Interestingly, the relaxedness parameter Xoff (bottom 
row of panels) behaves quite differently. It is relatively 
uncorrelated with MAH-PCi, but is fairly well corre- 
lated with MAH-PC2, at least for the full sample (since 
the relaxed sample is cut on Xofr, this makes it hard to 
tell how much correlation is present in this case). For 
the full halo sample, unrelaxed systems have a systemat- 
ically higher value of Xoff, and systems with larger values 
of Xoff have experienced significantly more acceleration 
in their MAH. 

5. APPLICATION TO OBSERVATIONS 
5.1. What Do Structural Properties Tell Us? 

Clearly the formation history of dark matter halos, 
summarized by their MAH, relates to their final struc- 
ture at a statistical level. Observationally, structure can 
be measured directly, whereas history cannot be. Thus 
a practical application of our results is to infer an unob- 
servable quantity from observable ones. The first ques- 
tion is, given a set of measured structural properties, 
what to they tell us about the MAH? 

From the analysis in Section 4.2, the parameter most 
closely and simply related to formation history is the 
concentration. But which aspect of the formation history 
does concentration trace most sensitively? Fig. 7 shows 
the correlation between concentration and a number of 
the age indicators introduced previously. These are 2:^, 
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Fig. 5. — Elongation E of the halo versus elongation of the progenitor particle distribution at 2; = 0.1,0.5, 1.0,2.4 for relaxed (left) and 
unrelaxed (right) halos. 
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Fig. 6. — Structural parameters versus MAH-PCs for the full sample (left-hand panels) and relaxed subsample (right-hand panels). 



the redshift by which the main progenitor of a halo had 
a fraction x of the final halo mass at z = 0, for x = 
0.2, 0.5, 0.8, and (M/Mo)^, the fraction of its final mass a 
halo had reached by redshift z. They represent the values 
obtained by intersecting the MAH Ai{z) with a series of 
horizontal or vertical lines respectively. Of all the age 
indicators, Z0.2 shows the tightest correlation; generally 
the correlation with Zx is tighter than with {M/Mq)z. 
Thus concentration is most closely related to the early 
part of a halo's formation history, when the dense central 
core of the halo is established. 

We can generalize the specific values of x considered in 
Fig. 7 to a continuous range of values, and test for cor- 



relations with other structural parameters. Fig. 8 shows 
how the strength of the correlation between a given struc- 
tural parameter and Zx varies as a function of x. We use 
the Spearman rank coefficient as our measure of correla- 
tion. Concentration is most strongly correlated with ^0.2, 
although the decrease in correlation at x = 0.1 may be an 
artifact of our MAH selection and or limited resolution 
at early times. Elongation E and sphericity c/a are rel- 
atively insensitive to x, with correlation peaking slightly 
in the range x =0.3-0.5. Relaxedness Xoff shows a simi- 
lar pattern of anti-correlation, with a broad peak around 
X =0.4-0.6. Spin and triaxiality are more weakly anti- 
correlated, with the significance peaking at intermediate 
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Fig. 7. — Concentration versus age indicators for the full sample. 
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Fig. 8. — Strength of the correlation between structural proper- 
ties and the formation redshift Zx-, the redshift by which the halo 
had built up a fraction x of its final mass at 2; = 0, as a function 
of X. 

or high X. We also see that mass M is anti-correlated 
with particularly for small values of x. This conveys 
the fact that halos with lower values of Z0.2 will be more 
massive on average, as expected from the conventional 
picture of hierarchical structure formation. 

5.2. Splitting Halo Samples by Structural Properties 

Given the scatter in all the correlations seen so far 
between structure and formation history, the measure- 
ment of structural parameters in any single system will 
give only a rough indication of its past history. A more 
interesting possibility is to average over sets of objects. 
As mentioned in Section 1, structural properties such as 
shape and concentration have now been measured for 
hundreds of galaxy cluster halos, and may soon be mea- 
sured reliably for galaxy halos in forthcoming weak lens- 
ing surveys. In the case of the most massive clusters, it 
may even be possible to constrain the full 3-dimensional 
shape by combining lensing with other observations (e.g. 
Morandi et al. 2011). In Fig. 9, we show how, by split- 
ting a sample of halos into subsamples based on struc- 
tural properties at z = 0, we can construct subsamples 



with systematically different formation histories, as in- 
dicated by the distribution of Z0.5 values. In the four 
panels we consider splits based on concentration, elon- 
gation, sphericity and relaxedness. Consistent with our 
earlier results, we find that a split based on concentra- 
tion produces the largest offset in the mean value of Z0.5 
between the two subsamples, while splits based on shape 
produce similar but smaller offsets. (These offsets will 
also be further reduced in most practical tests, since for 
most systems we will observe only projected shape, not 
true 3-D shape.) Relaxedness (which could be deter- 
mined from detailed mapping of the mass distribution, 
e.g. in the X-ray) does not produce much shift in zo.5 if 
we split the sample evenly, but it does select out system- 
atically younger halos if we cut out the minority of very 
unrelaxed systems. The cut illustrated in Fig. 9 is close 
to the one we have use throughout the paper to separate 
relaxed and unrelaxed systems. A similar selection by 
relaxedness has already been applied a pilot sample of 
10 clusters by Smith & Taylor (2008) in order to study 
the connection between X-ray morphology, offsets from 
X-ray scaling relations, and formation history, but larger 
samples are needed to reach definitive conclusions. 



Splits based on present-day samples can also tell us 
other things about the past history of a population of 
halos. In Fig. 10 we show how cuts based on elongation 
or sphericity at z = select halos formed from mate- 
rial with a systematically different spatial distribution 
at z = 0.5. This might eventually provide interesting 
tests of non-Gaussianity in structure formation, or test 
the dependence of galaxy formation on the detailed en- 
vironment around a group or cluster. 



6. SUMMARY 

The recent studies of Sll and Jll have clarified how 
the structural properties of halos are inter-related, and 
how they relate to a few specific measures of the for- 
mation history. Here we have extended this approach 
to study the connection between structure and forma- 
tion history in more detail. We first examine the intrin- 
sic variation in halo formation histories, as summarized 
by their MAHs A4{z). Applying principal component 
analysis to the individual steps of the MAHs themselves, 
we can decompose the range of formation history into 
an mean MAH and a main set of variations about the 
mean. The mean MAH and variations away from it along 
the first principal axis MAH-PCi are well-fit by the 2- 
parameter function suggested by McBride et al. (2009), 
although the two parameters f3 and 7 are not a natural 
parameterization of this variation. Since MAH-PCi is 
the strongest principal component, and since the vari- 
ations along this axis correspond to MAHs that reach 
a given fraction of the final mass at systematically ear- 
lier or later times, it is natural to consider this the best 
(MAH-based) definition of halo age. 

A second component, MAH-PC2, is also relatively 
important in the sense that it accounts for 17% of the to- 
tal variance in the MAHs. This component corresponds 
to acceleration or deceleration in the accretion rate at 
late times. It is not particularly well fit by the McBride 
formula. Third and higher PCs correspond roughly to 
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Fig. 9. — Distributions of the formation redshift 2:0.5 in halo samples split by present-day structural properties 
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terms in a Fourier-like decomposition of the MAH, with 
an additional oscillation about the mean in each succes- 
sive PC, but we have not studied them in detail since 
they are less important sources of scatter individually. 
Collectively, however, they account for almost 30% of 
the total variance, showing that individual MAHs are 
not particularly well fit by smooth curves. This is pre- 
sumably due to the large stochastic jumps produced by 
major mergers. 

Relating formation history and structural parameters, 
we recover the trends seen in the previous studies. Con- 
centration, in particular, correlates strongly with age 
indicators. Studying the relationship between concen- 
tration and MAH-PCi, however, we find that it is bi- 
modal. Relaxed systems have a fairly tight correlation 
between age (as expressed by MAH-PCi) and concen- 
tration, but in disturbed, low-concentration systems this 
correlation breaks down. Thus concentration provides 
an effective measure of age, but only in relaxed systems. 
Testing the strength of the correlation against various 
age-related parameters, we find it is tightest for the for- 
mation redshift 2:0.2, the redshift by which a halo had 
built up 20% of its final mass in its most massive pro- 
genitor, or equivalently the redshift where A4{z) = 0.2. 
Thus concentration is an indication of early merger his- 
tory. 

Testing for correlations with other parameters, we find 
that shape and relaxedness are better tracers of the more 
recent formation history of a halo. Shape parameters 
such as elongation or sphericity are reasonably well corre- 
lated with the shape of the mass distribution from which 



the halo formed, but only going back to low redshifts, 
that is to say only considering a halo's recent past. Simi- 
larly, irregularities in the density profile seem to correlate 
more with the last half of a halo's growth history. 

These theoretical results have several possible applica- 
tions to real systems. In the case of individual systems, 
concentration can give us a broad indication of early his- 
tory, while relaxedness can give us an indication of later 
history. Perhaps more interesting, however, are the im- 
plications for large samples of halos. Splitting cluster cat- 
alogs by concentration, shape, or relaxation should select 
sub-samples with systematically different formation his- 
tories. Comparing these to predicted MAH distributions 
could offer new cosmological tests, updating an early idea 
of testing cosmology by measuring the ages of clusters, 
as indicated by substructure (Richstone et al. 1992). On 
galaxy scales, if weak lensing surveys can reach the pre- 
cision where halo shape is measured routinely, this may 
allow a direct empirical test of how the MAH influences 
galaxy formation. 
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APPENDIX 

A. CORRELATION TABLES FOR HALO PROPERTIES 

Correlation plots and Spearman coefficients using the subset of 7 parameters also used by Jll is shown in Fig. Al, 
to aid comparison with this previous work. The full correlation analysis of the halo structural properties is shown in 



Fig. A2. 
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Fig. Al. — Spearman correlation table for the 7 parameters also used in Jeeson-Daniel et al. (2011). 
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Fig. A2. — A full Spearman correlation table of the parameters discussed in Section 2.6. The values indicated on the axes are the base-10 
log of the parameter in question. 



